FERMILAB-PUB- 10-028-E-PPD 



February 24, 2010 
Version 1.0 

Jet Physics at the Tevatron 

Anwar A Bhatti 
The Rockefeller University, 1230 York Ave, New York NY 10065 

Don Lincoln 

Fermi National Accelerator Laboratory, P.O. Box 500, Batavia IL 60510 

Abstract 

Jets have been used to verify the theory of quantum chromodynamics (QCD), measure the 
structure of the proton and to search for the physics beyond the Standard Model. In this article, 
we review the current status of jet physics at the Tevatron, a yfs = 1.96 TeV pp collider at the 
Fermi National Accelerator Laboratory. We report on recent measurements of the inclusive jet 
production cross section and the results of searches for physics beyond the Standard Model using 
jets. Dijet production measurements are also reported. 



I. INTRODUCTION 

The theory of quantum chromodynamics (QCD) [l| is currently the best description of 
the fundamental strong force. This theory describes the color interaction between quarks 
as being mediated by gluons, which are the vector bosons of the strong force. It has been 
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successfully tested in collisions between e + e |2j, ep [3|], pp |4j, and pp [5]. 

One of the basic properties of QCD is that its coupling strength a s decreases with the 
energy of the interaction and that, at sufficiently high energies, QCD calculations can be 
performed using perturbation theory in powers of a s |6[. Currently, these perturbartive QCD 
(pQCD) calculations are available at a next-to-leading order (NLO) for many processes and, 
in some cases, at next-to-next-to-leading order (NNLO) approximation. Leading order (LO) 
calculations, supplemented with parton shower calculations Ml, are used in several Monte 
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Carlo event generators |8|, |9[. In addition, matrix element generators which match NLO 
calculations of rates for QCD processes with a parton shower Monte Carlo event generator 
are also available [loj], but only a limited number of processes have been implemented. 

The pQCD calculations result in a small number of partons in the final state, while 
experimenters observe "jets" of particles. These jets retain the kinematic properties (energy 
and momentum) of the parent partons (quarks or gluons). In order to facilitate comparison 
between data and calculation, jet finding algorithms have been devised that are insensitive 
to the difficult-to-calculate low energy phenomena that govern the transition from low- 
multiplicity partons to high-multiplicity particle final states. There are several jet finding 
algorithms and the details of the measurements are sensitive to that choice. 

While events in which jets are created are used for a detailed understanding of the strong 
force, it is also possible that such events could also reveal new physical phenomena, including 
quark substructure (compositeness), extra spatial dimensions and new particles which decay 
into jets. Because of their high energy, jets can probe very small distances. At the Tevatron, 
the highest pt jets can probe distances down to O(10~ 17 ) cm. 

In 2001, the Fermilab Tevatron pp collider commenced its Run II, with a collision energy 
of 1.96 TeV. This energy is higher than the 1992 — 1996 Run I energy of 1.8 TeV. Even this 
relatively small increase in energy leads to a substantial increase in jet production with large 
transverse momentum, px, by about a factor of three at pt = 500 GeV/c. The beam intensity 
is much higher than Run I due to the addition of the Main Injector and the Recycler Ring 
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to the Fermilab accelerator complex. In addition, both the CDF [ll| and D0 detectors 
were upgraded. The results reported here utilize an order of magnitude higher integrated 
luminosity than reported previously ij. 



II. PERTURB ATIVE QCD 



The theory of QCD describes the behavior of those particles (quarks q and gluons g) that 
experience the strong force. It is broadly modeled on the theory of Quantum Electrodynam- 
ics (QED), which describes the interactions between electrically-charged particles. However, 
unlike the electrically-neutral photon of QED, the gluons, the force-mediating bosons of the 
strong interaction, carry the strong charge. This fact greatly increases the complexity in 
calculating the behavior of matter undergoing interactions via the strong force. 

The mathematical techniques required to make these calculations can be found in text- 
books (e.g. 13j). Instead of giving an exhaustive description of those techniques here, we 
focus on those aspects of the calculations employed most frequently in the experimental 
analysis, thereby clarifying the phenomena experimentalists investigate. 
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FIG. 1: Stylized hadron-hadron collision, with relevant features labeled. Note that a LO calculation 
of the hard scatter (dashed line) will assign a jet to final state radiation that would be included in 
the hard scatter calculation by a NLO calculation (dotted line). 
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FIG. 2: (a) Leading order Feynman diagrams, (b) Next to leading order loop diagrams, (c) Next 
to leading order tree diagrams. 

At high energies, the strong interactions between two hadrons can be factorized into 
three components: (a) the probability of finding the partons in the hadrons, (b) the in- 
teraction between quasi-free partons, and (c) fragmentation and hadronization of the final 
state partons. The process is schematically shown in Figure [TJ The cross section of the 
hadron-hadron scattering with 4-momenta Pi and P 2 can be written as 13[ : 

a(Px,P 2 ) = dxi dx 2 fi(x 1 ,fil)f j (x2,fi 2 F )a ij (p 1 ,p2,a s (fi 2 B ),Q 2 /fi 2 F ,Q 2 /fi 2 R ) (1) 
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The hard interaction between partons i,j is given by &ij where p\ = X\P\ and p 2 = x 2 P 2 
are 4-momenta of the two partons. It is independent of the incoming hadrons' structure 
and can be calculated using pQCD. The function fi(x,[i 2 F ) is the probability of finding a 
parton i with momentum fraction x at the scale /ip, and is called the parton distribution 
function PDF. The sum i, j is over the partons in the two respective hadrons. a s is the strong 
coupling constant. fiR is the renormalization scale, which is used to remove non-physical 



infinities inherent in fixed-order calculations. Q is the characteristic scale of the interaction 
and is related to some physical scale in the interaction, such as the of the leading jet. 

This factorization of hadron-hadron interactions into a short distance interaction (hard 
interaction) and a large distance interaction (PDFs) is done at an arbitrary energy scale 
/xp. A parton emitted below the scale /xp is considered to be part of the hadron structure 
and thus is described by the PDFs. The hadron structure is measured by many different 
experiments, especially lepton-hadron scattering experiments. The PDFs are determined 
by fitting these data and parametrized at a scale /Xq. The QCD evolution equations 
currently available at NNLO in perturbation theory, are used to evolve these PDFs to any 
different scale /x. The cross section is a convolution of the PDFs and the parton level cross 
section and, to be consistent, both must be calculated to the same order in perturbation 
theory. A complete calculation, including all orders in the perturbation series, must be 
independent of /xp, but a fixed order calculation depends on this choice. To evaluate the 
sensitivity to /xp, it is usually varied up and down by a factor of two. Scale dependencies 
are expected to decline with the addition of higher order terms in the calculation. 

In Figure EJ for jet production, examples of leading order (LO) and next-to-leader order 
(NLO) Feynman diagrams are shown. At NLO, the jet cross section receives contributions 
from virtual corrections to the two parton final state, and from real corrections from the 
three parton final state. Both contributions are divergent, but the sum is finite. One 
sees that a NLO pQCD calculation can describe up to three jets in the final state. For 
inclusive jet production, many different implementations of NLO pQCD calculations are 
available 14l-ll7l|. All these programs use Monte Carlo integration techniques to calculate 
the real and virtual contribution to the cross sections. The EKS program calculates 
the cross sections for user-predefined cuts on transverse momentum and rapidity ranges. 
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15] and NLOJet++[16| generate parton events with weights (both positive and 



negative) and thus full event kinematics are available to the user for jet clustering, detector 
acceptance calculations and study of any other distribution. These programs require a huge 
amount of CPU to reach the desired precision. A lar ge f raction of this CPU is consumed 
in evaluating the PDFs. In the FastNLO program [l7|, the convolution in Equation [1] 
is modified to a product of a perturbatively calculable piece ^ n ,i,k,i,m{^)i calculated using 
NLOJet++, the strong coupling constant a"(/x^ m ^), and a function F^Xa , %l\ /x^) which 
depends on only the PDFs and the factorization scale /x. The time consuming piece a is 
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calculated only once. The function F{ is calculated on a grid of different values of x± } X2, and 
fiF only once and interpolated between those points while evaluating the full cross section. 
This procedures significantly speeds up the calculations. 

Near threshold, the phase space for the emission of real gluons is limited and large 
logarithmic corrections to the above cross section calculation may arise from the incomplete 
cancellation of infrared divergences against the virtual gluon emission contributions. For 
jet production, these corrections are expected to contribute at very high x, where parton 
distributions are falling very steeply. For the inclusive jet cross section, these threshold 



corrections have been calculated to NNLO at next-to-leading-logarithmic accuracy 18 



and 



are found to be small. However, the corrected cross section shows a substantial reduction 
of the scale [1f,Hr dependence. 

The partons radiate when they pass through a color field. In this type of radiation, two 
forms dominate, collinear and soft radiation. Collinear radiation is in the direction of parent 
parton, while soft radiation is just low energy emission. These soft and collinear radiation 
can be calculated in a leading-logarithmic approximation to all orders and this algorithm is 
a crucial component of event generators. In this formulism, a gluon radiates another gluon 
or converts into a qq pair according to DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) 
splitting functions [7J which depend on a s and a variable z which describes the energy 
sharing between two daughter partons. Similarly, a quark or an anti-quark can radiate a 
gluon which radiates further. This sequential radiation results in a shower of partons. The 
radiation process is continued until the parton virtual mass t is smaller than a mass scale to- 
This simple procedure is augmented by angular ordering, i.e. each subsequent emission is 
required to have a smaller angle, to simulate color coherence effects. Color coherence leads 
to suppression of soft gluon radiation in certain regions of phase space. In the final state 
showers, the radiation is limited to a cone defined by color flow lines and the emission angle 
at each branch point is smaller than the previous emission angle. Partons with a virtual 
mass t < to ~ 1 GeV/c 2 are combined into hadrons using a phenomeno logical hadronization 
model. The hadronization models have been tuned to reproduce the jet structure observed 
at e + e~ colliders. The showering process described above is used for radiation from outgoing 
partons and called final state radiation (FSR). For radiation from incoming partons, initial 
state radiation (ISR), the event generators use so-called backward evolution. First the 
momentum fractions X\, 22 of partons participating in the hard interaction are determined. 



6 



Then, the parton shower that preceded the hard interaction is subsequently reconstructed, 
evolving partons from the hard interaction scale Q backward in time towards smaller Q 
where the PDF / is evaluated. The color coherence in initial state radiation is slightly more 
complicated but still follows the same angular ordering. Interference between initial state 
and final state radiation is implemented in herwig but not in pythia. The independent 
variable, t, which is how the evolution of a s (t) is parametrized, is not unique. In the pythia 
showering algorithm, the squared mass m 2 of the branching parton is used as the evolution 
variable, herwig uses t = m 2 /(2z(l — z)) where z = Ej,/E a is ratio of daughter parton 
energy (Ef,) to the parent parton energy (E a ). In recent versions of PYTHIA, the option of 
t = of the branching parton is also available as the evolution variable. 

Occasionally, the radiated parton is at a sufficiently large angle to the parent parton and 
carries enough energy that it leads to an identifiable jet. Because these jets are typically 
of lower £>T; it does not dominate the event kinematics. However, this lower radiation 
becomes important in studies of jet multiplicity. This part of ISR/FSR can be considered as 
a part of the hard scatter cr^- or treated independently. For instance, we see in Figured] two 
boxes surrounding the hard scatter. The dashed box surrounds the leading order scatter, 
while the dotted box surrounds a next-to-leading order diagram. At NLO pQCD, one parton 
emitted from either incoming or outgoing parton is part of the short distance hard cross 
section a^. In LO event generators, these parton emissions are treated quasi-independently 
of the hard scatter and are part of leading-log showering process. 

The partons from the incoming hadrons which do not participate in the hardest scatter 
in an event also interact, but these interactions are normally soft. The particles produced 
in these multi-parton interactions are, on average, isotropically distributed in the allowed 
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191 ] space and can overlap with the jets produced in the hard interactions. Being 



soft, multi-parton interactions are in the non-perturbative regime and thus are implemented 
in event generators using phenomenological models. The parameters of these multi-parton 
interaction models have been tuned to reproduce the transverse energy and multiplicity 
distributions of the particles observed far away from the hard jets in collider data 20, 



21] . The beam remnants, the partons which do not participate either in hard or multi- 



parton interactions, go along the beam direction. However, they do carry color and to 
become color-singlets, they must exchange (soft) partons with the rest of the event. The 
particles produced in multi-parton interactions and from the hadronization of beam remnants 
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collectively constitute the underlying event. 

As described above, pQCD predictions for jet production are available at NLO at the 
parton level only. These predictions can not, in principle, be compared directly with the 
data which is available at the particle level, because these parton level calculations do not 
include hadronization effects and the contribution from the underlying event. On the other 
hand, event generators include hadronization, underlying event energy, and ISR/FSR to 
all orders in the leading-log approximation, but the hard interaction is calculated at LO 
only. Thus to compare data with QCD predictions, a hybrid scheme is generally used. The 
parton level NLO calculations are correcting for the underlying event and hadronization 
effects before they are compared with data. These corrections are determined using Monte 
Carlo event generators by comparing the jets obtained by clustering the pQCD partons with 
the jets obtained by clustering the partons after the showering process. The parton jets are 
obtained from the Monte Carlo events in which the underlying event simulation (multiple 
parton scattering) has been turned off. 

Over the last decade, there has been a lot of progress in simulating high jet multiplicity 
events using tree-level matrix elements. In alpgen 22], events with the exclusive parton 
multiplicities n — 2, 3, 4, and 5 are generated using matrix elements from pQCD at the tree 
level. To include the effect of soft and collinear emission to all orders (albeit in the leading- 
log approximation), these events are passed through a showering program e.g. in pythia. 
The phase spaces of the matrix elements and the parton showering program overlap. In 
particular, showering programs occasionally generate hard partons which can lead to a state 
which has already been generated by the matrix element. To avoid this double counting, 
a matching criteria is used. For example, in alpgen it is required that the number of jets 
produced by clustering the partons produced by the matrix element is the same as those 
produced after the showering for n < 4. The events which do not satisfy this condition 
are rejected. For the n = 5 parton state, the showering algorithm is allowed to produce a 
higher jet multiplicity state. The spectra from each different multiplicity are combined to 
form the full spectrum. In sherpa 23], a different matching procedure 24] is used where 
parton showers above a cut off k^-like measure (c.f. Equation H]) are vetoed. Both alpgen 
and Sherpa have been extensively tested at the Tevatron in W/Z+jet production, but these 



studies are not discussed here due to space constraints 



25]. 



While lepton-nucleon deep inelastic scattering (DIS) experiments are able to precisely 
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measure the quark content of the proton, this precision is not achievable for the gluon, 
especially at high x. At low x, the gluon distribution can be determined precisely using QCD 
scaling violation in DIS data. Studying the high-x gluon distribution functions requires data 
from hadron-hadron scattering. The effect of including Tevatron jet data in global fits to 
determine PDFs is described in Section TV CI 



III. JET CLUSTERING ALGORITHMS 



Because a parton carries the strong charge, it is not directly observable. It showers into 
many partons which combine together to form a large number of particles which travel in 
roughly the same direction as the initial parton. The kinematic properties of the initial 
parton can be inferred either from the shower of partons or from the jet of collimated 
particles. For this inference, these particles or partons must be clustered into a jet by an 
algorithm. In pQCD, at NLO and higher orders, a jet algorithm is needed to define physics 
observables which are well-defined i.e. they are soft /collinear safe. Jet algorithms are run on 
a few partons generated in pQCD calculations to construct such variables. Experimentally, 
the final state particles are observed as tracks in the tracker systems or as towers of energy 
in the calorimeter. These tracks or towers must also combined into a jet so that they can 
be compared to the parton produced in the hard interaction. In the following, we will 
collectively call the (a) partons in a pQCD calculation, (b) partons or (c) particles produced 
in Monte Carlo event generators, or (d) towers or (e) tracks, or (f) reconstructed particles 
observed in a detector as the objects which are input for a clustering algorithm. 

For a valid comparison between observations and theoretical predictions, the clustering 



algorithm must satisfy some basic criteria 
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The algorithm must be safe against soft 



(infra-red IR) and collinear radiation, invariant to boosts along the beam direction, and 
should be insensitive to the non-perturbative hadronization effects. In an algorithm which 
is not safe against soft/collinear radiation, the virtual and real contributions in pQCD cal- 
culations do not cancel completely and thus the predicted cross sections are ill-defined. One 
should be able to run the same algorithm on the detector calorimeter towers or tracks, 
particles or the multi-parton state from event generators and partons in fixed order pQCD 
calculations and get sensible results. Experimentally, jet clustering algorithms should be 
insensitive to the energy from additional hadron-hadron collisions in the same bunch cross- 
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ing which overlaps the energy from hard interactions, and should not consume too much 
computer resources such as CPU. Finally, the algorithm must be completely specified to 
avoid different interpretations. 

Commonly used jet-finding algorithms can be divided into two categories: (a) cone clus- 
tering and (b) pair-wise recombination algorithms. With a few exceptions, only cone clus- 
tering algorithms have been used at hadron colliders. The cone clustering algorithms used 



prior to the Tevatron Run II were not IR/collinear safe 26|, |27|, and it was proposed to 
add an additional seed at the midpoint of stable cones. This made the new algorithm 
IR/collinear safe to NLO for the inclusive jet cross section measurement. For other physics 
observable, it is either safe at LO only or unsafe at all orders [281. Various issues related to 
jet reconstruction are extensively discussed in a recent review [29|. 

A. Cone Clustering Algorithm 

For jet studies in Run II, both the CDF and D0 collaborations are using the Midpoint 



algorithm, as laid out by the QCD Workshop recommendations [26[, but the two implemen- 
tations differ in some details. Below we describe the implementation of this algorithm by 
the CDF collaboration. 

The clustering process starts by making a list of all objects to be clustered. In simulated 
events, all the particles or partons are included without a p^ threshold. However, in data, 
the calorimeter towers are required to have pi > 100 MeV/c to minimize the effect of noise. 
From this list, a second list of seed objects is made with the requirement that the of the 
objects exceeds a fixed threshold of 1.0 GeV/c. At each seed location, the 4-momentum of 
the cluster is determined by summing the 4-momenta of all the objects within a distance 
R = (y — y r ) 2 + ((t> — (t> r ) 2 from the seed (y Cl (j) c ). The 4-momenta are summed using the 



E-scheme 26] 



(E, p x , p y , p z ) = ^2(E, p x , p y , p z )i (2) 

i 

Pr=v^ ^ = iln(§±£) ^ = tan-WpJ. (3) 



This scheme is different from the Snowmass scheme [30( used in Run I, where the clustering 
centroid was defined as the -fir-weighted average of r\ and <j). Using the center of the cluster 
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as a new seed location, the process is iterated until the center of the circle (y c , C ) coincides 
with the position of cluster 4-momentum. 

After all the stable cones have been identified beginning with real seeds, there is an 
additional search for stable cones using as seed locations the midpoints between the initial 
set of stable cones. The cone finding algorithm allows that the same object may be part 
of many cones. The shared objects are uniquely assigned to a single cone using a split- 



merge algorithm as specified in 



261 ] . If two stable cones share objects, the shared transverse 



energy is compared to the transverse energy of the lower px cone. If the ratio of the shared 
transveres energy to the transverse energy of the lower p^ cone is higher than the energy 
fraction j ' merge , the two cones are merged. Otherwise, based on proximity, the shared objects 
are assigned to the nearest cone. The two collaborations use different values of j merge- CDF 
(D0) uses f merge = 0.75 (0.50). This split-merge procedure may lead to jets which are not 
circular in y — <p space. After all the objects above threshold have been uniquely assigned 
to a stable cone, the jet kinematics are determined using the same E-scheme. 

A cone clustering algorithm can be made infrared safe to all orders if a stable cone is 
evaluated at each point in y — <p space. Such an algorithm is very CPU-intensive even when 
the number of particles is modest and thus is not practical beyond some parton level pQCD 
calculations. Recently, a new seedless cone clustering algorithm has been proposed which is 
infra-red and collinear safe to all orders in perturbation theory. The Seedless Infrared Safe 
Cone SISCONE algorithm 31] uses the fact that a circle enclosing a set of particles can be 
moved around such that two of the particles lie on its circumference. Consequently, all stable 
circles can be reconstructed by considering all possible pairs of particles. After determining 
all the stable circles h, the algorithm merges and splits the stable circles to uniquely assign 
the particles to a single circle. This algorithm is fast and has been used at the Tevatron for 
comparison purposes only. 



B. Pairwise Clustering Algorithm 



The cone algorithm combines all the objects within a distance R from the seed. In 
contrast, the recombination algorithms combine pairs of objects based on some measure dij 



and is an attempt to "undo" the showering of partons. The fcx algorithm 
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33| starts 



with a list of proto-jets given by 4-momentum (E,p x ,p y ,p z ). All the objects which are to 
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be clustered are considered as proto-jets. The transverse momentum p?, rapidity y, and 
azimuthal angle of a proto-jet are calculated using Equation [31 

For each proto-jet % and the pair i 7^ j), di and dij are defined as 

* = Ph <k = min(p%, P 2 l) te-ytf + ^-M 2 (4) 

where D is the parameter which controls the size of the jet. For the algorithm, the 
parameter p — 1. The algorithm determines the minimum d m i n of the di and all the dy. If 
dmin = di, the proto-jet is not mergable and is promoted to a jet. Otherwise, the proto-jets 
i,j are merged into a single proto-jet with the 4-momentum (Eij,pij) = (Ei + Ej,Pi + pj). 
The process is repeated until no proto-jets are left. 

The &t algorithm has been extensively used at e + e~ and ep colliders. At hadron colliders, 
the environment is more challenging. The energy from multi-parton interactions and beam 
remnants and pile-up can contribute to the jets and must be taken into account. The 
large particle multiplicity observed in hadron-hadron collisions requires substantial CPU 
resources to process an event. Thus the use of the k^ algorithm has been limited at hadron 
colliders. The D0 collaboration measured the inclusive jet cross section in Run I 34j. In 
the Tevatron Run II, the fcx algorithm has only been used by the CDF Collaboration to 
measure the inclusive jet cross section (described in Section IV Bl) . 



Recently, two more recombination algorithms using p = (Cambridge- Aachen) 35| and 



p — — 1 (anti-A/p) 36| in Equation H] have been proposed. The Cambridge-Aachen algorithm 
combines particles based only on their relative distance. The anti-^x algorithm combines 
the highest p^ objects in the events first. This leads to circular jets, which have well-defined 
area like the cone jets. Thus far, these algorithms have not been used at the Tevatron. 



IV. JET ENERGY SCALE DETERMINATION 



At the Tevatron, jets are generally measured using a calorimeter, which is sensitive to 
both charged and neutral particles. Both CDF and D0 utilize sampling calorimeters, which 
measure only a small fraction of the energy of the particles. This observed energy is multi- 
plied by a calibration constant so that it is equal to the sum of the energies of the incident 
particles. The calorimeter response is different for hadrons, photons, electron and muons. 
For hadrons, the response depends on the momentum and the flavor of the particles, whereas 
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for photons and electrons it is almost momentum-independent. Muons normally deposit a 
little energy (~ 1 GeV) in the calorimeter, which is almost independent of the muon mo- 
mentum. Neutrinos escape without interacting and lead to an imbalance in the measured 
Px hi the event. The observed jet energy must be corrected for the calorimeter response 
and other detector effects. The two collaborations employ different techniques to determine 



these jet energy scale corrections. The CDF collaboration's technique [37J depends on an 
accurate modeling of the calorimeter response to single particles and a knowledge of the 
Px spectrum of the particles in a jet, whereas the D0 technique is data-driven and utilizes 
the fact that in photon-jet events pr,jet = Pt,, [38[. These techniques are used to calibrate 
the central region of the calorimeter where the tracking system is available to measure the 
charged particle momentum and also the calorimeter response is uniform. 

This approach was applied in the optimum calorimeter region for both collaborations. 
The calorimeter response was extended to other regions (0.1 < \r]\, \r]\ > 0.7 for CDF and 
\rj\ > 0.5 for D0), by using dijet balancing to scale the jet energy response in the other 
regions to the one in the optimum region. The energy from additional pp interactions in the 
same bunch crossing is subtracted, based on the number of reconstructed primary vertices 
in an event. For cone jets, this correction is determined from minimum bias events by 
summing the energy in towers in a cone of radius R placed randomly in the calorimeter. 
The procedure for fc T jets is described in Section IV Bl 

a. Photon-Jet balancing: In this technique, the jet energy is determined by scaling 
the measured jet px to the photon's px in photon-jet events. The photon energy is measured 
by the electromagnetic (EM) calorimeter which is linear and has very good energy resolution. 
In the approximation of 2 — > 2 scattering, the jet transverse energy px is equal to the photon 
Pt- The real situation is a little more complicated due to presence of initial state radiation 
ISR, the energy not clustered in the jet, and contributions to the clustered jet from the 
multiple parton interactions. To be insensitive to these effects (especially ISR), D0 evaluates 
the missing px ($r) projection fraction along the photon direction using: 

R had = 1 + (5) 

The hadronic recoil correction factor, Rhad, is the scale factor to the entire recoil system. By 
requiring that the jet is back-to-back with the photon and, in the absence of any additional 
jet(s) in the event, Rhad is almost equal to the jet response. The derived response is expressed 

13 



in the jet energy E' determined from the pt of the photon and the position of the balancing 
jet using E' = p\j sin ^ et as both p\ and the direction of the jet are accurately measured and 
thus E' provides a better estimate of the jet energy than the direct jet energy measurement 
by the calorimeter. It is preferred over jet pp> as the calorimeter response depends on the 
energy of the incident particles and thus parametrization of calorimeter response in E' is 
more natural. 

The EM calorimeter is calibrated using the electrons from Z boson decays such that the 
reconstructed Z boson mass is equal to the world average 39( . The EM calorimeter response 
to electrons and photons is similar, but not the same, as photons start their shower later 
than electrons. This difference is small and is evaluated using simulated events at pr = 100 
GeV/c. The estimated uncertainty on the photon energy scale is 0.5% at low E' and 0.8% 
at high E'. Using this procedure, the D0 collaboration has achieved a 1% accuracy on the 
jet energy scale in photon-jet events. The current statistics of the 7+jets sample limit the 
direct measurement of the jet energy corrections in the central region to E' < 350 GeV. The 
response is extrapolated to higher energies using Monte Carlo, which has been tuned to the 
data. The correction to a single jet with a given algorithm and size is deduced from Rhad 
using simulated events. 

The calorimeter response to jets depends on their flavor, as the particle spectrum and 
multiplicity for quark-initiated and gluon-initiated jets are different. The jet energy scale 
corrections determined from 7+jet events is valid only for the flavor composition of 7+jet 
events. Event topologies with different flavor composition will have different jet energy scale 
corrections. Indeed, D0 tuned the single-pion response in their detector simulator to data 
and used pythia to generate photon-quark, photon-gluon and dijet events. They found 
that the gluon jet response was 8(2)% lower for jets with 20(500) GeV of energy. In QCD 
jet production, the fraction of gluon-initiated jets changes with jet p^ and the corrections 
were adjusted to account for this variation in flavor composition. With these additional 
corrections, the uncertainty on the jet energy scale is reduced to an unprecedentedly-small 
value. 

b. Jet Corrections using Single Particle Response: Another approach to deter- 
mine the jet energy correction is based on a knowledge of the calorimeter response to each 
particle that makes up a jet. The CDF collaboration measured the calorimeter response to 
charged hadrons and electrons using both pp collider and test beam data. The calorimeter 
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simulation was tuned to reproduce the measured response. The calorimeter response to a 
jet was determined by a convolution of the single particle response with the type and mo- 
mentum distribution of particles constituting a jet as given by a fragmentation model. CDF 
used QCD dijet events with the pythia fragmentation model to measure the default jet 
corrections. The pythia fragmentation model agreed well with the particle pp and multi- 
plicity distributions in a jet measured in pp data. The herwig event generator was used to 
crosscheck the pythia fragmentation functions and the results determined using two gener- 
ators were found to agree well. In this procedure, the difference in calorimeter response to 
gluon-initiated and quark-initiated jets is automatically included. Although this procedure 
requires a detailed knowledge of the calorimeter response and a well-tuned simulation, it has 
the advantages that the correction can be easily determined for any event topology over the 
entire kinematic range, and real and simulated data have the same corrections and thus can 
be treated on an equal footing. 



V. INCLUSIVE JET CROSS SECTION 



The inclusive jet cross section measurement [40h45|] has been used to test QCD and to 
search for physics beyond the Standard Model by searching for an excess of events at large pp. 
During Run I at ^/s = 1.8 TeV, the search was limited by both theoretical and experimental 
systematic uncertainties even with an integrated luminosity of 100 pb _1 . The uncertaint y on 
the jet energy scale dominated the experimental uncertainty. NLO pQCD calculations 
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161 ] significantly reduced the dependence on the factorization and renormalization scales and 
the remaining dependence is ~ 10% and almost independent of jet pr for pr > 100 GeV/c. 
The jet cross section is not very sensitive to non-perturbative hadronization effects. The 
underlying event contributes approximately 2 GeV/c of pp to a jet and thus is significant 
only at low pp. The main theoretical uncertainty arises from uncertainty in the parton 
distribution functions, especially for large x gluons. Unfortunately, inclusive jet production 
is the only process in which the high x gluon distribution can be directly measured. The 
other possible process, photon-jet production, has a limited reach in x and the associated 
theoretical uncertainties are large. The gluon distributions are also measured from lepton- 
proton scattering data through QCD scaling violations, but these measurements are also 
limited to low x values. Because of these limitations, the Run II inclusive jet cross section 
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has been primarily used to constrain the gluon content of the proton. The data at high y are 
particularly useful as it probes high x at lower Q values where the contribution of physics 
beyond the Standard Model, if any, is negligible. 

Due to the higher center of mass energy and much larger integrated luminosity, Run II jet 
measurements extend the jet spectrum to higher jet pp compared to the Run I measurements, 
by approximately 200 GeV/c. Both collaborations implemented an improved jet clustering 
algorithm. The jet clustering algorithms used in Run II are IR/collinear safe at least to 
the order of the available pQCD calculations. An accurate determination of the calorimeter 
response from the pp data and also a refinement in the techniques to determine the jet 
energy scale have lead to reduced uncertainty compared to Run I. In previous inclusive 
measurements, the hadronization effects were ignored as they were much smaller than both 
the experimental and theoretical uncertainties. In Run II, both collaborations have evaluated 
the effect of hadronization and corrected the parton level pQCD calculation. In addition, 
the pQCD calculations are corrected for the energy from the underlying event, determined 
using tuned event generators. In contrast, in Run I energy from the underlying event was 
removed from jets in data. 

The CDF collaboration measured the inclusive jet cross section using the cone cluster- 
ing 46] with cone size R = 0.7 and the clustering algorithm with D = 0.5, 0.7 and 
1.0 47]. The D0 collaboration has recently published the inclusive jet cross section 
using cone clustering algorithm with cone size R = 0.7. These three measurements are 
described below. 



A. Measurement using Cone Clustering Algorithm 

1. D0 Collaboration 

The D0 collaboration analyzed 0.7 fb" 1 of data taken during 2004 — 2005 to measure the 
inclusive jet cross section for px > 50 GeV/c in six rapidity bins, |Ay| = 0.4 wide, over the 
range < \y\ < 2.4. The data were collected by triggering on a jet passing a px threshold. 
Six triggers with p^ thresholds of 15, 25, 45, 65, 95, and 125 GeV/c were used to collect data. 
Due to high production rates, only predetermined fractions of the lower threshold triggers 
were recorded. The efficiency for triggering on jet events was measured using data collected 
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FIG. 3: (top) D0's observed inclusive jet differential cross sections corrected to the particle level 



48(. The NLO 



in six rapidity regions compared to next to leading order (NLO) QCD predictions 
QCD predictions are calculated with the CTEQ6.5M parton distribution functions; (bottom) ratios 
of the measured cross sections over the NLO QCD predictions. The data agrees with the theory 
quite well and has remarkably small systematic uncertainties. 



with a muon trigger which did not rely on calorimeter activity. These different jet triggers 
were combined to form the full Pt spectrum with each trigger contributing to a unique Pt 
range. Only those data for which the trigger efficiency is > 98% are used. The events were 
required to have a reconstructed primary vertex and the position of the pp interaction be 
within 50 cm of the detector center along the beam direction. This requirement ensured 
that the jets follow the projective geometry of the calorimeter and thus their energy was 
accurately measured. The consequence of this requirement was a reduction of only 7.0±0.5% 
in the integrated luminosity. The primary vertex was reconstructed using charged particle 
tracks measured using silicon micro-strip and scintillating fiber detectors located inside a 
solenoidal magnetic field of 2 T [ljj]. 
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The triggered data includes events containing cosmic ray interactions, beam halo and 
detector noise. These contributions are mostly asymmetric and lead to a large imbalance 
in the momentum in the plane transverse to beam direction, ftp. In contrast, for QCD jet 
production, is ideally zero, apart from a small neutrino contribution. In QCD events, frp 
arises mainly from fluctuations in calorimeter response and is much smaller than the total 
energy observed in the detector. Most of these background events are removed by requiring 
the ratio of fir to the transverse momentum of the leading jet to be small. Remaining 
backgrounds are removed by requiring that the shape of energy deposition in the calorimeters 
be consistent with the expected shape from a hadronic jet. The shape of energy deposition 
for a jet is very different from the energy deposited by a cosmic muon or a beam halo 
particle, as a jet consists of many particles. These shape requirements also remove photons 
and electrons. These requirements are highly efficient for the signal and the remaining 
background is estimated to be < 0.1%. 

The measured j>t of each jet is corrected for calorimeter non-linearity and energy lost 
in uninstrumented regions. These average jet-by-jet energy corrections do not correct the 
smearing (bin-to-bin migration) of jets due to the finite energy resolution. This smearing is 
determined using an iterative procedure. It is assumed that the particle level physics (true) 
spectrum is described by the function 

JPf \ AT ( P T Y fl 2 P T COsh (ymin) \ P , , , * 

F{P ^ y) = N0 {T00GeVr c ) 7i J eXp( ~ 7PT) ' (6) 

where y m i n is the rapidity lower bin edge. This functional form is a good representation 
of the NLO pQCD prediction and fits the measured raw inclusive jet spectrum well. This 
true spectrum is smeared using the jet energy resolution function, which is determined us- 
ing pp collider dijet data and simulated dijet events. The resulting smeared spectrum is 
compared with data using a \ 2 test- The process is iterated to determine the best param- 
eters (N ,a,{3, 7) of the true function, F(p T ,y). This true spectrum is used to correct the 
migrations between bins in px in the observed data. In the central region, the migration 
correction is a multiplicative factor that is 0.8 — 0.9 at low p^ and 0.7 at higher p^, with 
a strong dependence on y. The true spectrum F(pr r ,y) is measured separately for each ra- 
pidity bin. The jet rapidity is measured very precisely and thus migration between rapidity 
bins is small. The y migration corrections are less than 2% in most bins and 10% in the 
highest p T bin where spectrum is the steepest. The rapidity unsmearing is applied after the 
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Pt unsmearing. After the jet energy scale and resolution smearing corrections, the observed 
data distribution has been corrected to the particle-level jets and is completely independent 
of the detector properties. 

The observed inclusive particle jet spectrum compared with the NLO pQCD predictions, 
corrected for underlying event and hadronization effects, is shown in Figure [3J The NLO 



pQCD predictions are calculated using NLO Jet++ and FastNLO 16|, LL7| which use ma- 
trix elements with /ir = fip = p 1 ^ . The parton distribution functions from CTEQ6.5M 49| 
are used, which include the Run I inclusive jet data. The dashed curve shows the NLO pQCD 
prediction calculation using MRST 2004 divided by the same calculated using CTEQ6.5M. 
The experimental uncertainty, dominated by the jet energy scale uncertainty, is 12% at 
Pt = 50 GeV/c and 17% at p T = 550 GeV/c in the \y\ < 0.4 bin. The uncertainty is 
higher in other y bins. These results are the most precise to-date. The main theoretical 
uncertainties arise from the uncertainty on the PDFs and the missing higher order terms 
in the perturbation series. As is customary, the effect of higher order terms is evaluated by 
varying the renormalization and factorization scale. Fortunately, the change in cross section 
from varying these scales is almost independent of the pt of the jets. The predicted cross 
section changed by ~ 10 — 15% when the scale is changed to /i = 2p^ or fi = p^/2. There is 
a good agreement between the data and the theoretical predictions over the whole px range 
which spans 50 GeV/c to 550 GeV/c. Over this pt range, the cross sections falls by 10 
orders of magnitude. 

The data prefers the lower bound of the theoretical prediction, favoring a smaller gluon 
content of the proton at high x. The theoretical uncertainty arising from the uncertainties in 
the parton distribution functions is larger than the experimental uncertainties. These data, 
along with CDF inclusive jet data, have been used in the global fits to improve the precision 
of the gluon distribution function. The results of these fits are described in Section IV C[ 

2. CDF Collaboration 

The CDF collaboration measured the inclusive iet cross section for a cone size of R = 

I — I 

0.7, using slightly more data, corresponding to 1.13 fb _1 [461 ] . The measurement spans five 



rapidity bins \y\ < 0.1, 0.1 < \y\ < 0.7,0.7 < \y\ < 1.1, 1 



These bins are matched to the CDF calorimeter structure H( and thus are different than the 



< \y\ < 1.6, 1.6 < \y\ < 2.1. 
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binning used in D0 analysis. The data was collected between 2002 — 2005 using jet triggers 
with four thresholds: 20, 50, 70 and 100 GeV/c. In order to not saturate the data acquisition 
system by jet triggers, only 1/808, 1/35, 1/8 of lower threshold triggers were recorded. The 
transverse energy of each jet is corrected on average to form the jet p^ spectrum which 
was corrected for bin-to-bin migration of jets due to finite jet energy resolution. CDF 
used simulated events to evaluate the smearing corrections. The corrections depend on 
the shape of the true jet px spectrum and the jet energy resolution. A large sample of 
QCD jet events was generated using the pythia event generator §J and passed through 
the CDF detector simulation. The detector simulation was tuned to describe the single 



particle response measured in pp collisions [37] • These simulated data were analyzed using 
the same procedure as the one used for the real data to obtained the smeared spectrum. 
The bin-to-bin migration effect was determined by taking the ratio of smeared jet spectrum 
and the particle jet spectrum. For this procedure to be valid, the smeared spectrum of 
the simulated events must match the spectrum measured in data. The two spectra are very 
close but not exactly same. The simulated particle jet smeared pt spectrum was adjusted 
(re-weighted) to force it to agree with the measured spectrum. Re-weighting changes the 
unsmearing corrections by only a few percent. The unsmearing correction is < 5% for 
< 300 GeV/c and increases to as much as 20% at pr = 500 GeV/c. 
The corrected je t yy spectrum is compared to perturbative QCD predictions evaluated 
with the FastNLO 17J program using the CTEQ6.1M parton distribution functions 50 ]. 
The renormalization and factorization scales (hr,hf) are chosen to be pr/2, which are the 



same as used in the global QCD analysis to determine the PDFs 50|. Using — /If — Pt* 
gives up to 10% smaller predictions in the cross section. The perturbative QCD predictions 
are corrected for underlying event and hadronization effects measured using the procedure 
described in Section [Til While clustering the partons produced by the FastNLO, CDF used 
an ad-hoc parameter R sep which was introduced to mimic the split and merge procedure in 
iterative cone clustering algorithms 5 If . At order af, the final state can have up to three 
partons. Depending on their relative pt and their separation in y — space, these partons 
are clustered into two or three jets, Two partons are clustered into a single jet if they are 
within R from the jet centroid and within R x R sep of each other. A value of R sep = 1.3 is 
used in this calculation. An R sep = 2.0 (i.e. the midpoint algorithm without R sep ) yields less 
than a 5% increase in cross section for NLO QCD predictions. As shown in Figure IU the 
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data are in good agreement with the theoretical predictions. The experimental uncertainty, 
dominated by the jet energy scale uncertainty, is comparable to the theoretical uncertainty 
which is dominated by the PDF uncertainty. 
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FIG. 4: (left) CDF's observed inclusive jet differential cross sections corrected to the particle level 



46(. The NLO 



in five rapidity regions compared to next to leading order (NLO) QCD predictions 
QCD predictions are calculated with the CTEQ6.1M parton distribution functions, (right) Ratios 
of the measured cross sections over the NLO QCD predictions. The theory describes the data quite 
well. 



B. Measurement using the fcx Clustering Algorithm 

The clustering algorithm, which combines objects in pairs to reconstruct a jet, is infra- 
red and collinear safe at all orders in perturbation theory and is preferred over iterative cone 
algorithms. The CDF collaboration has measured j^] the inclusive jet cross section using 
the fop clustering algorithm with D = 0.4,0.7, and 1.0 using 1.0 fb -1 of data in the same 
rapidity region as used in the cone-based analysis 46j. The analysis procedure is similar to 
the one described in Section IV A2\ except the correction for multiple- interactions which are 
determined using a novel approach. The measured jet transverse momenta are corrected for 
this effect by removing a certain amount of transverse momentum, 8™ x (Ny — 1) where 
Ny denotes the number of primary vertices in the event. The value of 5™ was determined 
by requiring the shape of the pp spectrum at high instantaneous luminosity (£ mst ) to be the 
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same as the one at the low £ mst . After making the shape of two spectra the same, the data 
at low £ mst and high £ mst are combined. The study was carried out independently for each 
rapidity region and the results were consistent with a common value = 1.86±0.23 GeV/c. 
The corresponding correction for the cone jets is 0.97 ± 0.29 GeV/c, which is measured by 
summing the p^ in a cone of R = 0.7 in minimum bias events. 

The pt spectra are compared in Figure [5] with NLO QCD predictions using the 



CTEQ6.1M PDFs 



50] with /i = 0.5 x p™ axJct for D = 0.7. The theoretical predictions 



are calculated using the JETRAD [15[ program. The data are in very good agreement with 
QCD predictions except in the highest rapidity bin (1.6 < \y\ < 2.1), where the data are 
lower than the prediction, but well within experimental systematic uncertainties. The theo- 
retical uncertainties are dominated by the PDF uncertainties and are comparable or larger 
than the experimental uncertainties. The theoretical predictions using the MRST2004 PDFs 
is very close to those based on the CTEQ6.1M PDFs, except in the 1.6 < \y\ < 2.1 bin, where 
the MRST2004 cross section is smaller, but within the PDF uncertainty on the CTEQ6.1M 
prediction. The results for jet size D = 0.4 and D = 1.0 show similar behavior. 

The jet pt spectra measured using two different clustering algorithms are expected to be 
different and can be compared only via theoretical predictions. The ratios of data/theory 
from two analyses were compared and the two ratios were in very good agreement with each 
other except in the 0.7 < \y\ < 1.1 region where the /ct cross section is ~ 5% higher. In 
this y region, the CDF calorimeter coverage is not uniform which leads to a large variation 
in calorimeter response and poor jet energy resolution. The two CDF analyses have similar 
experimental uncertainties. Thus one concludes that both the ki and the cone clustering 
algorithms can be successfully used at the hadron colliders. 

C. Determination of Gluon Distribution Function 

The parton distribution function (PDF) fi(x, fi), which is the probability to find a parton 
with a type % = g,q,q with momentum fraction x and mass scale /x, must be experimen- 
tally determined. The PDFs for gluons and light quarks and anti-quarks (u, d, s, u, d, s) are 
normally determined from experimental data. For heavier quarks, i.e. c and b, they are 
normally dynamically generated through gluon splitting. Data from e ± p collisions at the 



ZEUS and HI experiments up, up, un, un collisions at CCFR/NuTeV [ 



and Drell-Yan 
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FIG. 5: (top) CDF's inclusive jet differential cross sections measured using the &t clustering 
algorithm in five rapidity regions. The data has been corrected to the particle jets and is compared 



471 ] . The NLO QCD predictions are calculated 



to next to leading order (NLO) QCD predictions 
with the CTEQ6.1M parton distribution functions, (bottom) Ratios of the measured cross sections 
over the NLO QCD predictions. The perturbative NLO QCD calculations are in good agreement 
with the data. 



(lepton pairs and W/Z bosons) production in pp 53] and pp collisions 54J, jet data from 
Tevatron and data from many other experiments, especially low energy, are used to extract 
the PDFs using a global fit. These experiments are sensitive to different fi. For example, 
e ± p experiments are sensitive to the sum of the q and q distributions weighted by e 2 q ^ and 
can not distinguish between quark and anti-quark distributions. Neutrino and anti-neutrino 
data are used to differentiate between q and q. The Tevatron jet data play a significant 
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role in constraining the gluon distribution at large x. The gluon distribution at low x are 
mainly determined from the scaling violations in lepton-nucleon scattering data. Normally, 
the results of these global analysis are fit at some initial scale Hq using tens of parameters. 
The PDF can be evolved to any arbitrary scale \i using QCD evolution equations which are 
available at NNLO approximation in perturbation theory. 



Both the MSTW 



551 ] and CTEQ 56[ collaborations have included Run II jet data in 



their global analysis to update the PDFs. The three Run II inclusive jet measurements are 
more accurate than Run I measurements, span a larger range and are consistent with 
each other 55]. As two CDF measurements used the same data, the MSTW collaboration 
decided to use the jet spectrum whereas the CTEQ collaboration is using the cone-based 
measurements. The Run I jet measurements do not play a significant role in the fit and 
thus the MSTW collaboration has dropped those data from the new fits. Comparisons of 
the gluon distribution g(x) determined in the new MSTW fit with the gluon distributions 
from MRST2004 53] and CTEQ6.6 58j] fits along with the 90% uncertainty band are shown 
in Figure M (left) for /i 2 = Q 2 = 10 4 GeV 2 . The new g(x) is lower than previous fits for 
x > 0.3 but within the still large systematic uncertainties. As a s and g(x) always appear as 
a product, the values of a s and g(x) are strongly correlated. The value of a s in the three sets 
of PDFs is different and thus g(x) is also expected to be slightly different. The fractional 
uncertainty on the gluon distribution is shown in Figure |6] (right). At x — 0.4 and /i 2 = 10 4 
GeV 2 , the uncertainty reduces from 18% when the jet data are excluded from the fit to 12% 
when jet data are included in the fit. This modest extra constraint will make the predictions 
more precise at the LHC in processes where gluon-quark scattering dominates. 



D. Determination of the strong coupling constant 

The hard cross section for jet production (aV, ) in Equation (JTJ at large pr can be expanded 
in powers of the strong coupling constant to the n th order in perturbation theory 

m=0 

where the perturbative coefficients are functions of the kinematic variables and the 
factorization and renormalization scales only. The coefficients are available for m = 
and m — 1. Using the jet cross section measurement, the strong coupling constant a s can 
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FIG. 6: (left) Ratios of various gluon PDFs to the MSTW2008 PDF. (right) The MSTW2008 



gluon PDF uncertainties for variations on the inclusion or exclusion of Tevatron data 



55]. 



be determined using Equation [T] provided the parton distributions functions are known. 
This technique was used by the CDF collaboration to measure a, at different px values 
and show its running with the hard scattering scale in Run I 59[. In Run II, the D0 



collaboration has used the same principle but an improved technique to measure a s [60] from 
the data used to measure the inclusive jet cross section. The hadronization and underlying 
event corrections were determined using pythia and applied to theory predictions. The 
hadronization (underlying event) corrections vary between -15% (+30%) to -3% (+6%), for 
p T = 50 GeV/c to 600 GeV/c [48]. 

The perturbative results are the sum of O(a^) pQCD calculation Q, supplemented 
with (9(o4)(2-loop) corrections for the threshold effects 18|. The PDFs are taken from the 
MSTW2008 next-to-next-to-leading order (NNLO) parametrization and the renormaliza- 
tion and factorization fiR,fiF scales are set equal to the pt of the jet. The NNLO (NLO) 
MSTW2008 PDFs are available for 21 different values of a s (M z ) ranging from 0.107 - 0.117 
(0.110 - 0.130) in steps of 0.001. This a s is used to evolve PDFs. 

Commonly available parton distributions from the MSTW and CTEQ groups include 
Tevatron jet data in the global fit. To avoid any correlation between the input PDFs and 
the extracted a s , only 22 of 110 available jet data points are used. These selected points 
contribute to the x region (x < 0.25) where PDFs in global fits are mainly determined by 
other experimental data and are not strongly influenced by Tevatron jet data. The jet data 



starts to affect g 
5% for x < 0.25 



x) at x ~ 0.2. The change in g(x) due to inclusion of jet data is less than 
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The central a s (Mz) result is obtained by minimizing \ 2 with respect to a s [Mz) and 
integrating over the nuisance parameters for the correlated uncertainties. The variation of 
a s (px) v s Pt is shown in Figure [7] (top). The running of a s as a function of jet follows the 
QCD evolution equation. The data points from the HI and ZEUS experiments follow the 
same curve but have large uncertainties. The a s , evolved to \x = Mz, is shown in Figure [7] 
(bottom). The combined result of 22 selected points a s [Mz) = 0.1161 



with the world average 0.1184 ± 0.0007 
are large. 
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FIG. 7: a s measurement as a function of the transverse momentum of the jet from several ex- 
periments, as compared to the expected variation of a s , setting a s (Mz) = 0.1184 ± 0.0007 [601 ] . 



VI. SEARCH FOR PHYSICS BEYOND THE STANDARD MODEL USING JETS 

Due to large theoretical and experimental uncertainties in the jet production rate, the 



high tail cannot be used to improve upon the current limits on new interactions [62|. 
However, there are kinematic distributions from jet events that can be used to identify 
physics beyond the Standard Model. Below we describe these searches using the dijet mass 
spectrum and the dijet angular distributions. These searches are not very sensitive to the 
jet energy scale, parton distribution functions or renormalization or factorization scales. 
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A. Dijet Mass Spectrum 



Many new physics models predict particles which decay into two high px jets. These 
particles can be identified by the reconstructed mass of the dijet system provided their 
intrinsic mass width is narrow. Such models include excited quarks 63], axigluons 



flavor-universal colorons 
heavy vector bosons 



65], color-octet techni-p [66], Randall Sundrum (RS 



68] and diquarks in the string-inspired E 6 model 69]. The excited 



gravitons 



6j, 



67j, 



quarks q* decay into qg. Heavy vector bosons W, Z' decay into qq or qq'. The axigluon 
A decays into qq, and E§ diquarks D (D c ) decay into qq(qq). The RS Graviton G* and 
color-octet techni-p p^ 8 both decay into either a qq or gg pair but their branching ratios are 
different. 

All these models predict an intrinsic mass width which is much smaller than both the 
detector resolution and mass broadening effects due to QCD radiation. These models can 
be divided into three categories depending on the decay channel, i.e. gg, gq and qq. The 
expected mass shapes for q* , G* , W, and Z' particles with a mass of 800 GeV/c 2 are 
shown in Figure [SJ Because q* and G* decay into gluons, their widths are broader than the 
widths for W and Z' . Gluons radiate more than quarks, resulting in a broader dijet mass 
distribution. These distributions are close and change the final limits by only 10-20%. These 
shapes can be used to search for resonance structure, independent of the model details. 
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FIG. 8: Expected dijet mass distributions for simulated signals for the following new physical 
models: q* — > qg, RS graviton (— > gg, qq) and W' — > qq and Z' — > qq with a mass of 800 GeV/c 2 . 



To measure the dijet mass spectrum [70j, the CDF collaboration used the same data set 
as used in the inclusive jet cross section measurement described in Section TV A 21 46] . The 
dijet mass is reconstructed from the two highest pt jets using 
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FIG. 9: (left) CDF's measurement of the dijet mass cross section for events in which the two 



highest pt have \y\ < 1.0 



701 ] compared to NLO calculations using CTEQ6.1M. (right) Ratio of 



data to NLO theory. The experimental uncertainties are dominated by jet energy scale uncertainty 
and are comparable to the theoretical PDF uncertainty. 



Jets produced by new physics are expected to be produced more centrally than by Standard 
Model processes and we expect a better signal to background ratio in the central region. 
Thus only those events in which two leading jets have \y\ < 1.0 are used. Moreover, the 
CDF calorimeter is best understood in this region. As shown in Figure [91 these data, after 
all corrections, are in good agreement with the NLO QCD calculations. The dijet mass 
spectrum before smearing corrections is shown in Figure [TO] Smearing corrections are not 
used to avoid any degradation of the resonant structure, if any, in the data. The measured 
mass spectrum is fit to a smooth background given by 

drrijj 

where po,Pi,P2, and are free parameters. The dijet mass spectra predicted by pythia, 
herwig, and NLO pQCD can be described well by this functional form. The fit to the 
measured dijet mass spectrum is shown in Figure fTOT a). The data are well described by this 
smooth function with a y 2 of 16 for 17 degrees of freedom. The deviation from the smooth 
curve is shown in Figure HOTb). These data are used to determine the exclusion limits on 
the existence of new particles decaying into jets, as there is no evidence for the existence 
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FIG. 10: (a) CDF's measured dijet mass spectrum [70J]. The dashed curve shows the fit to Equa- 
tion [9l Also shown are the predicted dijet mass distributions of the excited quark q* for the masses 
of 300, 500, 700, 900 and 1100 GeV/c 2 respectively, (b) The fractional difference between the 
measured dijet mass difference and the fit (points) compared to q* signals divided by the fit to the 
measured dijet mass spectrum (curves). 

of any resonant structure. The experimental limits are determined for the a s%9 = a ■ B ■ A, 
where a is the theoretical new particle production rate, B is probability of its decaying into 
two jets and A is the kinematic acceptance of the resulting particle jets to have \y\ < 1.0. 
The upper limits on a s%9 are evaluating using a likelihood function 

L = H^exp(- f i l )/n l \ (10) 

i 

where /ij = nf 9 +nf CD is the predicted number of events in bin i. The QCD dijet background 
nf D is determined using Equation Q by evaluating • • Arrijj ■ da/drrijj\i where Am 3 j, 
€i, and Ci, are the bin width, the trigger efficiency and the integrated luminosity for bin i 
respectively. The expected signal events nf 9 is given by a s%9 ■ Ci - ej • (ni/n to t) where is the 
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signal event selection efficiency in the i th dijet mass bin and ni/n to t is the predicted signal 
fraction in bin i. For each value of a sl9 , the likelihood is maximized with respect to the four 
parameters in Equation (jH]). This profiled likelihood is integrated over Bayesian priors for 



Til ] . A flat prior on a s%9 is used to 



the parameters describing the systematic uncertainties 
extract a Bayesian upper limit on that parameter. 

The mass exclusion limits for W', Z', q* and G* are determined by comparing the limits 
obtained using their respective signal shapes and the predicted theoretical cross section. For 
other models, the limits obtained for the above four signal shapes are used. The q* signal 
shape (qg) is used for axigluons, the flavor-universal coloron and the Eq diquark, as these 
particles do not decay into modes which include gluons and thus their signal shapes are 
expected to be narrower than the q* signal shape. For pxs, the limits obtained for the G* 
shape are used. Both the prs and the G* decay into qq or gg, but the branching fraction of G* 
into gg is higher. Thus in all the above cases, the obtained exclusion limits are conservative. 
The limits for these models are given in Table [B 

TABLE I: The limits on the masses of particles decaying into dijet in various models. Limits are 
in units of GeV/c 2 [70|. 
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The dijet mass spectrum measured by the D0 collaboration in six rapidity bins 72j is 
shown in Figure [TTJ The rapidity bin is labeled by the higher of the two jet rapidities. 
The data are compared to NLO pQCD predictions computed by FastNLO 17J using the 
MSTW2008 NLO PDFs [551 with scale PR = = (pt,i +p T ,2)/2 and are corrected for 
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non-perturbative effects determined using the pythia event generator. The data and QCD 
predictions are in reasonable agreement for the \y\ < 0.4 region which is not surprising as 
the MSTW2008 global fit includes the inclusive jet data (see Section IV CD . For higher \y\ 
bins, the data are below the theoretical predictions but within 1 a of the total experimental 
systematic uncertainty. The D0 collaboration is searching for new particles using these 
data. 



B. Dijet angular distributions 

The angle between the initial and final state partons in the center of momentum frame 
is sensitive to the spin of the exchanged or the intermediate particle and thus can be used 
to search for physics beyond the Standard Model. At hadron colliders, dijet production is 
dominated by the t channel exchange of a gluon, a massless vector boson, and the angular 
distribution has the familiar Rutherford scattering form 

d ° 1 = 1 (n) 

d cos 9* ~ (1 -cos9*) 2 sin 4 (fl*/2) 1 ' 

where 9* is the angle between the jet and the beam direction in the dijet center of momentum 
frame. The angular distribution of the new particles proposed in many new physics scenarios 
is relatively flat in cos 9*. For example, the angular distribution of spin 1 particles (W, Z', 
Axigluon, coloron) decaying in fermions is da jd cos 9* ~ l + cos 2 9*. Theories in which quarks 
are composite particles but with the compositeness scale much hi ghe r than the available 



energy, can be parametrized by an effective Lagrangian of the type [62j, 



g 2 

£ = Vj^KQil^Q^Vljl^W i = L,R, j = L,R 



where A is a parameter in the theory which controls the characteristic energy of the new 
interactions. The parameter t] is ±1 and determines the sign of interference between new 
interactions and the SM interactions. The main effect of substructure is to increase the 
proportion of centrally produced jets, which can be observed in the jet angular distribu- 
tions [73]. 



The ADD LED models Q,EJ proposed to solve the hierarchy problem, i.e. the difference 
between electroweak scale (~ 100 GeV) and the Plank scale Mpi (~ 10 19 GeV), assume 
the existence of extra spatial dimensions in which gravity is allowed to propagate. As a 
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consequence, gravity appears weak in the three conventional spatial dimensions. The Planck 
scale, the number of extra dimension n, their size R and an effective Planck scale Ms are 
related by Mpi = MsR n . Experimentally, Ms can be measured for different values of n. 
The Kaluza-Klein excitations of the graviton can be exchanged between partons and thus 
contribute to jet production, resulting in jets which are central. There are two different 
formalisms to describe LED models, GRW and HLZ 78|. In the HLZ formalism, the 
sub-leading dependence on the number n of extra dimensions is also included. 
In some models 
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81] . extra dimensions are assumed to exist at the TeV 1 distance 



scale. In these models, Kaluza-Klein excitations of SM bosons modify various production 
cross sections. In these models, gluons can propagate through the extra dimensions, which 
changes the jet cross section. The strength of the interaction is given by the model parameter, 
the compactification scale, Mc- 

To search for new physics, instead of studying the cos6>* distribution directly, it is con- 
venient to use the Xdijet distribution which removes the Rutherford singularity: Xdijet is 
defined as expQyx — 2/2 1 ) where y\ and 1/2 are rapidities [19[ of the two highest p<r jets in an 
event. For 2 — y 2 scattering of massless partons, the variable Xdijet is related to the partonic 
center-of-momentum frame polar angle 8* by Xdijet — (1 + cos#*)/(l — cosd*). 

The CDF collaboration studied the dijet angular distributions using 106 pb -1 of data 



from Run I [82j. The data excludes at 95% CL a model of quark substructure in which 
only up and down quarks are composite and the contact interaction scale is A^ d < 1.6 
TeV or A~ d < 1.4 TeV where the subscript refers to the flavor of quarks assumed to be 
composite and the superscript ± refers to the sign of the interference. For a model in which 
all quarks are composite, the excluded regions are A + < 1.8 TeV and A~ < 1.6 TeV. In 
Run II, the D0 collaboration measured the Xdijet distribution in the Xdijet < 16 range in 10 
dijet mass mjj bins covering the 0.25 < rtijj < 1.1 TeV/c 2 range using up to 0.7 fb -1 of 



data collected during 2004 - 2005 



83[ . The boost of the two-jet system is required to be 



Z/boost = 0.5 x I yi + y2 1 < 1. This requirement, combined with the dijet mass cut and the 
range of the Xdijet distribution, restricts the highest allowable rapidity to 1 2/1,2 j < 2.4 where 
the D0 detector performance is well understood. The measured distributions are corrected 
for detector effects using events generated with PYTHIA v6.419 js| with tune QW 21] and 
the MSTW2008LO parton distributions functions. This procedure corrects for the migration 
between dijet mass bins, as well as the shape of the Xdijet distributions in each mass bin. 
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The corrected normalized differential cross section distributions (1/crdijet • da/dxdijet) at the 
particle level are shown in Figure [12] for 10 dijet mass bins. The NLO pQCD predictions 
are computed using FastNLO [17j based on NLOJet++ [16j. These parton level predictions 
are corrected for the hadronization and underlying event contributions, which are evaluated 
using pythia. The theoretical uncertainties on the SM Xdijet distributions arising from the 
uncertainty on the PDFs and the uncertainty on renormalization and factorization scales are 
less than 2% and 5% respectively. These data are in good agreement with the SM predictions 
and thus are used to set exclusion limits in the parameter space of quark compositeness, 
ADD LED and TeV" 1 models. Calculations for all these models are available only at leading 
order, while pQCD calculations can be performed at next to leading order. For this analysis, 
the expected distributions for each new model are calculated at LO and then scaled by k- 
factors (k = ctnlo/clo) determined from pQCD calculations. The /c-factors vary from 1.25 
to 1.5. All these models predict a higher rate as Xdijet ~~ >" 1 an d as rrijj increases. However, 
the magnitude of the excess is different for different models. A Bayesian procedure 84J is 
used to obtain 95% CL. limits on the mass scale parameters A, Mc and Mg in the above 
models. The results in which the prior is chosen to be flat in the model cross section are 
given in Table [Til Other choices give similar but slightly higher limits 83 1 . 

The limits on Mc obtained in this analysis are the first direct search for TeV -1 extra 
dimensions at a particle collider, though inferior to indirect limits from precision electroweak 
measurements 81] . The limits on Ms in different formalisms of ADD LED are on average 
slightly higher than the recent D0 results obtained using dielectron and diphoton data 85 |. 
The quark compositeness limits are the most stringent limits to date. 



VII. CONCLUSIONS 



Since the start of the Run II, there has been a significant increase in the experimental 
data used in the jet analyses at the Tevatron. Due to the higher production cross section and 
the increase in available integrated luminosity, the inclusive jet cross section measurement 
has been extended to transverse momenta of 600 GeV/c. The experimental uncertainty 
is still dominated by the uncertainty on the jet energy scale which is +31/-26% in the 
Pt = 457 — 527 GeV/c bin in 0.1 < \y\ < 0.7 region for the CDF measurements and +16.0/- 
15.5% in the pr = 490 — 540 GeV/c bin in \y\ < 0.5 region for D0 measurements. Both 
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TABLE II: Expected and observed 95% C.L. limits on various new physics models. 
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FIG. 11: (left) D0's dijet mass measurement for six rapidity bins [72| compared to NLO using 
MSTW2008 PDFs. (right) Ratio of data/theory. The experimental systematic uncertainties are 
very small. The dashed line shows the effect of using the CTEQ6.6M PDFs. 
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FIG. 12: D0's measurement of the dijet x f° r t en dijet mass bins [83|. The solid line is the 
Standard Model prediction at NLO. The dashed lines show the predictions for various new physics 
models. 



collaborations used a modern jet finding algorithm, the midpoint algorithm, which is infra- 
red and collinear safe at next-to-leading order in pQCD for measurements of the inclusive jet 
cross section. In addition, the treatment of the underlying event energy and hadronization 
effects has improved over the techniques used in Run I. The experimental uncertainties are 
lower than the theoretical uncertainties, which are dominated by the uncertainties on the 
parton distribution functions. These data have been used in global fits to determine the 
parton distributions and have decreased the uncertainty on the gluon distribution for the 
x > 0.3 region. The new gluon distribution is slightly lower than that determined by Run I 
jet cross section measurements. This increase in the accuracy of gluon distributions will 
make the prediction of various processes more precise. 

These jet data have been successfully used to search for physics beyond the Standard 
Model using jet kinematic distributions. The dijet mass spectrum has been used to expand 
the exclusion regions in parameters of excited quarks in quark compositeness models, E 6 
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diquarks, axigluons and heavy vector bosons W and Z bosons, and the techi-p in color 
octet models. The 95% CL. lower mass limits range from 630 GeV/c 2 for colorons to 1.25 
TeV/c 2 for axigluons. The dijet angular distribution has been used to extend the limits 
on the quark compositeness mass scale, ADD large extra dimensions and the TeV -1 extra 
dimensions. The 95% CL. lower limit on the compositeness mass scale is 2.8 TeV. The 
95% CL. lower limit for compactification mass scale in the TeV -1 model is 1.5 TeV. The 
limits on the ADD large extra dimensions range from 1.9 TeV to 1.3 TeV depending on the 
number of extra dimensions in the HLZ formalism. The 95% CL. in the GRW formalism is 
1.6 TeV. In most cases, these are the best limits to date. 

We thank Andrew Beretvas, Kenichi Hatakeyama and Marek Zielinski for reading the 
manuscript and for insightful comments. 
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